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Abstract 

i-C It is shown that the classical dynamics of 1:1 resonance interaction 

q between two identical linearly coupled Duffing oscillators is equivalent 

| ; to the symmetric (non-biased) case of 'macroscopic' quantum dynamics 

of two weakly coupled Bose-Einstein condensates. The analogy develops 
through the boson Josephson junction equations, however, reduced to a 
single conservative energy partition (EP) oscillator. The derived oscilla- 
tor is solvable in quadratures, furthermore it admits asymptotic solution 
£h in terms of elementary functions after transition to the action-angle vari- 

ables. Energy partition and coherency indexes are introduced to provide a 
i ^ | complete characterization of the system dynamic states through the state 

variables of the EP oscillator. In particular, nonlinear normal and local 
r- H mode dynamics of the original system associate with equilibrium points 

^ of such oscillator. Additional equilibrium points - the local modes - may 

occur on high energy level as a result of the symmetry breaking bifurca- 
tion, which is equivalent to the macroscopic quantum self-trapping effect 
in boson Josephson junction. Finally, since the Hamiltonian of EP oscil- 
lator is always quadratic with respect its linear momentum, the idea of 
second quantization can be explored without usual transition to the rigid 
pendulum approximation. 

<N 

1 INTRODUCTION 

The present study brings attention to the analogy between 1:1 resonance dynam- 



ics of coupled classical Duffing oscillators and macroscopic quantum dynamics 
of coupled Bose-Einstein condensates PJ. Despite of quite different physical 
contents the resultant equations for both types of interaction appear to have 
surprisingly the same form of boson Josephson junction (BJJ) tunneling equa- 
tions [27] [25] . These represent a classical single degree-of- freedom Hamilto- 
nian system whose conjugate states are the fractional population disbalance 
and phase difference. In macroscopic quantum dynamics, such equations are 
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eventually obtained from the nonlinear Schrodinger (Gross-Pitaevskii) equation 
[5] by representing the corresponding wave function as a superposition of 
two wave functions [27] [25]. In the case of classical resonance dynamics, the 
BJJ equations are obtained after specific coordinate transformation by apply- 
ing parameter variation and conventional averaging techniques to the differential 
equations of motion of interacting oscillators. Further analogies are revealed by 
establishing the links between qualitative features of the BJJ equations and 
the corresponding physical effects developed in both classical and macroscopic 
quantum dynamics. In particular, specifics of the energy exchange between non- 
linear oscillators due to their 1:1 resonance and so-called nonlinear normal mode 
motions are discussed first. It is known that, during the normal mode motions 
[2"B] [3"U] , the energy partition between the oscillators from one cycle to another 
is fixed so that there in no energy flow between the oscillators. However, when 
the initial states of oscillators are out of compliance with any of the normal 
modes, the oscillators slowly exchange by some portion or even all of the energy 
in a beat-wise manner. Such beating phenomena have been in focus of nonlinear 
physics and physical mechanics for few decades by very different theoretical and 
practical reasons [5] [§]. Usually, the beat dynamics are described through the 
transition to amplitude-angle or similar coordinates followed by averaging the 
new equations over one cycle of vibration. The resultant system may appear 
to have one or may be several integrals dictated by the system' symmetry [B]. 
These types of integrals actually provide the description for the energy exchange 
between the interacting oscillators. Recently, based on the detailed parametric 
study of phase trajectories, the notion of limiting phase trajectory was intro- 
duced [13] [14] [12]. It was noticed that, in the limit when the entire energy of 
the system swings from one oscillator to another, the descriptive phase angles 
resemble state variables of impact oscillators with non-smooth temporal shapes. 
The importance of this observation is that it provides asymptotic simplifications 
for the cases extremely opposite to the normal modes. Further, analytical algo- 
rithms of nonsmooth temporal transformations [20] [21] where adapted to build 
approximate analytical solutions for limiting phase trajectories [T7] [TB] [5S]. In 
the present work, it is shown that the phase variable, which determines the en- 
ergy partition (EP) between the two oscillators, is given by a strongly nonlinear 
conservative oscillator, which is exactly solvable for any intensity of the energy 
exchange. Furthermore, the EP oscillator admits two different asymptotic limits 
associated with the vanishing or maximal possible intensities of beats in which 
the oscillator becomes harmonic or vibro-impact, respectively. In particular, the 
impact limit corresponds to the limiting phase trajectory on which the energy 
exchange involves the total energy of the system. Finally, it is shown that first 
integral of the EP oscillator completely determines necessary and sufficient con- 
ditions for the nonlinear mode localization phenomenon [26] |30| . It is shown 
that vibrations near nonlinear local modes of coupled oscillators are equivalent 
to the macroscopic quantum self-trapping (MQST) effect. As discussed earlier 
[2"5] . there are two types of MQST, such as running-phase MQST and 7r-phasc 
MQST. Note that the BJJ captures both types of MQST whereas the super- 
conductor Josephson junction (SJJ) rigid pendulum approximation ignores the 
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7r-phase MQST. An early overview and introduction to SJJ can be found in 
PP. According to this theory, a single wave function associates with a macro- 
scopic number of electrons which are assumed to condense in the same quantum 
(superconductive) state. When two superconductors are close enough to each 
other, practically at about 30 A, the superconductors begin to interact through 
the tunneling effect such that at about 10 A between the superconductors their 
macroscopic quantum phases cannot be viewed as independent any more and 
must be described as a single system pp. In a similar way, it is impossible 
to separately consider two even weakly coupled classical oscillators when they 
resonate. This is due to the fact that slowly oscillating resonance energy flows 
between the oscillators cannot be ignored anymore and must be interpreted as 
a new system. As mentioned above, such new system represents a strongly 
nonlinear conservative oscillator with quite interesting physical and mathemat- 
ical properties. Surprisingly enough, the basic case of such oscillator has been 
known for several decades as unique example of strongly nonlinear exactly solv- 
able oscillator with no relation to any physically meaningful situation [7] , except 
few phenomenological applications [TH] [3] [22] ■ Sections [2] through [5] describe 
the methodology based on the classical model of coupled Duffing oscillators. 
In Section [6] it is shown that the resultant 1:1 classic resonance equations take 
the form of BJJ equations. Then, the BJJ derivation in macroscopic quantum 
dynamics is briefly recalled for comparison, reason. Finally, the analogies and 
physical meaning of the BJJ variables in classical and quantum dynamical cases 
are discussed. 

2 NONLINEAR BEATS 

2.1 Coupled nonlinear oscillators 

Let us consider a system of two identical linearly coupled unit-mass oscillators 
(Fig.l) described by Hamiltonian 

H = \(vi + vl) + n(m) + n> 2 ) + ~& ( Ul - u 2 f (i) 

where Ui and Vi (i = 1, 2) are the coordinates and linear momenta, respectively, 
H(ui) is an even analytic function describing the potential energy for each of the 
two oscillators, and b is the coupling stiffness, which is assumed to be relatively 
week, 6/n"(0) < 1. 

Introducing the parameters, Q — y/H"(0) + b and e — bQ~ 2 , and subtracting 
the parabolic component from the potential energy, 

U(u) - ^n"(0)ii 2 = eU(u) (2) 

brings Hamiltonian ([I]) to the form, which incorporates the additional assump- 
tion that both nonlinearity and coupling are of the same order of magnitude, 

H = X - (v\ + vl) + ^il 2 (u\ + u 2 2 ) + e [U( Ul ) - n 2 Ul u 2 + U(u 2 )] (3) 
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Note also that the coupling is represented now in somewhat canonical form 
after the non-coupled terms of the interaction energy Q, buf/2, have been asso- 
ciated with the corresponding oscillators. The differential equations of motion 
are given by 

dH . dH 

Ui = —, v t = -— (4) 

OVi OUi 

or 

Ul = Vl 

u 2 = v 2 

in = -n 2 ui+e[n 2 u 2 -U'{ui)} (5) 
v 2 = -n 2 u 2 + e[n 2 u 1 -U / (u 2 )} 

As e — > 0, system ^ degenerates into two identical harmonic oscillators 
whose total energies are separately conserved. At non-zero e, the oscillators 
become non-linear and coupled with each other in such way that one of the 
oscillators is loaded proportionally to the displacement of another oscillator. 
Since system ([5]) is perfectly symmetric and conservative, it is reasonable to 
assume a relatively slow energy exchange between the oscillators due to the 
weak coupling. In order to describe the energy exchange dynamics in phys- 
ically meaningful terms, let us introduce a new set of variables as follows 
{«!,«!, «2,«a}- > {K(t),6(t),5(t),A(t)}: 

n? f e A 

ui = v K cos — I — cosd 

Ui = — \fKVlcos\ — | — | sin 8 
\2 4/ 

u 2 = -^Bta(| + £W + A) (6) 
v 2 = VKQsin ( - + - ] sin((5 + A) 



4, 

If K, 8, and A are constant, and S = fit, expressions ^ represent the exact 
general solution of the decoupled set of harmonic oscillators (j5J), e = 0. There- 
fore, relationships ^ implement the idea of parameter variations compensating 
the perturbation, when £ ^ 0. In order to track the oscillator energies during 
the vibration process, let us introduce quantities 

Ei = ^(v 2 + n 2 u 2 ) = ^E o (l-sm0) 

E 2 = l -(v 2 + n 2 u 2 )= l -E a {l + sine) (7) 
Ei 2 = -{v\V2 + tt 2 uiu 2 ) = -^JE 1 E 2 cos A 



where 



E a = Ei+E 2 = ^n 2 K (8) 
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Expressions |7]) and ([8| clarify physical meaning of the variables K, 9 and 
A participating in transformation ([6]), while the variable 5 is the fast phase 
associated with the principal temporal rate of the vibrating oscillators. In par- 
ticular, the quantity K is proportional to the total energy of the decoupled and 
linearized oscillators, the phase 9 characterizes the energy partition between the 
oscillators, and A relates to the phase shift between the oscillators as discussed 
below. In case £ ^ 0, the energy parameter K has small temporal fluctuations 
due to the coupling and nonlinear terms in ([5]). Nevertheless, expressions and 
([8| can still be used for characterization of the energy distribution between the 
oscillators. For that purpose, as follows from (|7|), the interval — 7r/2 < 9 < n/2, 
on which Eq > E\ > and < E 2 < Eq, is sufficient. In particular, the case 
9 = corresponds to equipartition, E\ = E 2 , under which the system oscillate 
either out-of-phase (A = 0) or in-phase (A = tt) according to the sign conven- 
tion in ([6]). It is convenient to deal with one full period of the phase shift A 
within the range — ir/2 < A < 3tt/2, including both in-phase and out-of-phase 
modes. 



2.2 Energy partition and coherency indexes 



As an alternative to the angular quantities 9 and A, let us introduce the energy 
partition index P and coherency index Q describing the system' vibrating states 
as follows 

Ei~E 2 . n (9) 



P = 



E 1 + E 2 



= — sm t 



- 1< P < 1 



Q = 



E 



12 



= -cosA, - 1 < Q < 1 



(10) 



v 7 E~iE~2 

According to definition ([9| , the number P = indicates the energy equipar- 
tition, Ei = E 2 , whereas P = 1 or P = - 1 correspond to the case when all 
the energy belongs to either first or second oscillator, respectively. According 
to definition (10 1, the boundaries Q = 1 and Q = — 1 correspond to in-phase 
and out-of-phase modes, respectively, as illustrated by Figs. 1 and 2. 

The derivation of equations for the descriptive phase variables includes sub- 
stitution of the coordinate transformation (|6| in |5]), which is in fact the varia- 
tion of constants procedure, and then implementing the averaging with respect 
to the fast phase S as described in Appendix. Assuming the monomial form, 
U{ui) = auf/A, such averaging gives 



K = 
9 = en sin A 

— eJ7(cos A tan 9 — n sin ( 



A 

5 



(11) 



= n 



1 



cos A tan 



k{1 - sin0) 



where the following nonlinearity parameter is introduced 

_ 3aK 



(12) 
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The first equation in (11) shows that the energy parameter K remains av- 



eragely constant regardless the perturbation parameter e. The fact that K is 
constant justifies the use of quantities ^ and ^ for characterization of the en- 
ergy exchange between the oscillators since neither the coupling nor nonlinear 
stiffness in ^ can accumulate the energy during one vibration cycle. In or- 
der to clarify physical meaning of the parameter k, consider a single oscillator, 
ii + Q 2 u + au 3 = 0, whose mean (over the period) potential energy components, 
corresponding to linear and nonlinear stiffness terms, are Eq = il 2 < u 2 > /2 
and E a — a < u 4 > /4, respectively. Assuming the harmonic temporal mode 
for the coordinate u(t), and taking into account ^ and (12), gives 



Eg 

En 



(13) 



Therefore, k characterizes the strength of nonlinearity in terms of its relative 



energy capacity during one vibration cycle. Taking into account (11) and (12), 
gives also k — 0. As follows from (11), further complete description of the 



dynamics can be conducted now in terms of the two phase shift parameters, 
A(i) and 9(t); the fast phase 8(t) is obtained then by integration from the last 
equation in pT] ). 

As discussed above, the corresponding dynamical system on the phase plane 
A — 9 will be considered on the rectangular 



D = {-it/2 < A < 3tt/2, -tt/2 < 9 < tt/2} 



(14) 



In the low nonlinearity range, < K < 1, there are two stationary points, 
(0, 0) and (tt, 0), corresponding to the out-of-phase and in- phase vibration modes, 
respectively. However, two more stationary points (± arccos K , 0) emerge from 
zero when the nonlinearity becomes sufficiently strong, 1 < k; see Fig. 1. Geo- 
metrical interpretation of the bifurcation at k = 1 will be given below in Section 

o 



3 EP OSCILLATOR 



It can be shown by inspection that system (11) admits the integral as follows 

(15) 



G = — cos A cos ( 



1 



- k cos 9 — const. 



Taking into account ( 15 ) and eliminating the phase A from the second and 



third equations of system (|11|, gives a single strongly nonlinear conservative 

(16) 



oscillator with respect to the coordinate 9 in the form 
9 + (efl) 2 [G 



n f ,n tan 9 1 o „ , 

2 1 ^ 2 - — k sin 19 = 



cos 2 9 



Note that the quantity G in equation ( 16 1 remains constant only on a fixed 



dynamic trajectory in the plane #-A but may vary from one trajectory to an- 
other. Therefore, the number G ( 15 ) must be calculated first by fixing some 
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point {#o,A } on the trajectory. Then equation (16 1 can be solved by making- 
sure that the initial condition {0(0), 0(0)} corresponds to the fixed trajectory 
according to the second equation in ( 11 ). The parameter k however can be cho- 
sen regardless the dynamics in the plane 9- A. According to ([7]), equation (16 1 
constitutes a principal equation describing the energy exchange between oscilla- 
tors ([5]). Moreover, if the function 9(t) is known, then other two phase variables, 
A and 6, are obtained from system (|11[) by differentiation and integration. 



It is shown therefore that the entire system (111 is exactly solvable in quadra- 



tures since general solution of equation ( 16 ) can be obtained from its 'energy' 
integral Hg = const, 



H e = ~e~ + {eQ) z ( ^G 2 tan 



1 i» 

2 



1 

16 



k 2 cos 29 



(17) 



The dynamics of oscillator ( 16 1 essentially depends on the shape of its poten- 



tial energy within the interval —tt/2 < 9 < tt/2. In particular, if the parameter 
k is small enough, then oscillator (17 1 has one stable equilibrium position at 
9 = 0. However, high energy levels of the original system 0, that is large k, 
can make the equilibrium position 9 = unstable by generating two new (sta- 
ble) equilibrium positions. Such kind of bifurcation obviously relates to that 
described at the end of Section [2j According to definition |9|, every equilibrium 
of the oscillator ( 16 1 corresponds to a fixed (no beats) energy partition between 



the oscillators in the original system, in other words, - nonlinear normal mode. 



Recall that oscillator ( 16 1 is obtained after the procedure of averaging has been 
applied to the original model For validation purposes, Fig. 3 illustrates 

the behavior of energy partition index ^ based on numerical solutions of both 
types of models the original model ^ and oscillator (16). The initial configu- 
rations are chosen to cover the areas of normal and local modes as illustrated 
by the model trajectories on configuration planes in Fig. 4. Overall, Figs. 4 
through 7 illustrate the dynamic trajectories on configuration planes and the cor- 
responding behaviors of energy partition ^ and coherency ( 10 ) indexes, when 
increasing the energy (nonlinearity) level above the critical point k = K* = 1. 
Such diagrams give quite complete characterization of the dynamics by showing 
how the energy is distributed between the oscillators and what is the phase shift 



between the oscillators at any given time. Note that the derived oscillator ( 16 1 



captures both local dynamics near separate nonlinear normal modes and drifts 
over multiple modes, including stable and unstable ones. This essentially com- 
plements the theory of nonlinear normal modes in classical dynamics covering 
mainly periodic normal mode motions and their small neighborhoods |30j . 
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4 ASYMPTOTICS OF BEATS 



4.1 Linearized case 



When the original system is linear (a = k = 0), equation (161 admits 
explicit analytical solution within the class of elementary functions 



0(t) 



arcsm sin 



?o sin 



(18) 



where 9q is the amplitude of 9, whereas another constant can be introduced 
into the phase <f>(t) = e sec o \G\£lt as an arbitrary temporal shift admitted by 
equation ( |16[ ). As mentioned in Introduction, such type of explicit solution 
has been known for quite a long time with no relation to any physical system 
[7], however it was used recently in some physical and mechanical applications 
[19], [3] in a phenomenological way. Although solution (181 holds only for the 
linear model, k = , it nevertheless helps to clarify specifics of the behavior of 



phase variables in nonlinear cases. In particular, substituting ( 18 ) in the second 



equation of (11), gives 



A(i) 



arcsm 



\G\ tan 6*o cos <fi(t) 
y 1 - sin 2 6»o sin 2 <j>(t) 



(19) 



Fig. 8 illustrates the relationship between the beat dynamics of the linearized 
system and the corresponding phase variables, 9(t) and A(t). The coordinates 
Ui(t), (i = 1,2) represent exact analytical solution under the initial conditions 
obtained from Q at t = 0. The integral G is calculated at the amplitude poinlQ 
9q = 7r/2— 0.01, chosen slightly below its maximal possible magnitude ir/2. Since 
A = as 9 = 6q, then G — — cos#o and therefore </>(£) = eVLt. Other parameters 
are taken as e — 0.01, SI = 1.0, K = 1.0 and 5(0) = 0. As follows from Fig. 
8, the behavior of phase variables 9 and A resembles smoothed time histories 
of the coordinate and velocity of a simple impact oscillator. As mentioned 
in Introduction, this fact was noticed first in [T2], [T3] based on the analysis 
of phase equations similar to (111 however obtained in a different way after 
complexification of the original coordinates. In particular, it was found that 
the 'impact limit' corresponds to the most intensive energy exchange between 
the oscillators when each of the oscillators periodically hosts the total energy of 
the system. It is seen now that, in the linearized case, such asymptotic follows 



directly from exact solutions, ( 18 ) and ( 19 ) 



9(t) 
A(t) 



. . , n f'2 

arcsm(sm <p) = —t 



arcsm 



cos < 



| cos < 
eflt, 9 -> tt/2 



7T 



7T 

= 2 6 



(20) 



lr The notation 6o should not be confused with the initial value, which is 9(0) = 0, according 
to the present form of the solution. 



8 



where t(z) and e(z) are triangular sine and rectangular cosine wave functions 
whose amplitude is unity and the period is normalized to four in order to provide 
the basic relationships of nonsmooth temporal transformations |21j . t'(z) = e{z) 
and e 2 (z) = 1. 

Now substituting (18) in Q, gives 



Ei = -Eq(1 — sin O sineQt) 



En 



-E (l + sin 6*o sinefit) 



(21) 



Taking into account ^ and |2l| ), gives the corresponding energy partition 
index 

P{t) = - sin0 o sineOt (22) 

Recall that the number P = indicates equipartition, Ei — E%, whereas 
P = 1 or P = —1 correspond to the case when all the energy belongs to either 
first or second oscillator, respectively. As follows from (22), such states can be 
reached only in the limit case (20), when 8q =%/2 . 



In case of small amplitudes \9 \ -C n/2, corresponding to a moderate energy 
exchange, solutions ( 18 ) and ( 19 ) are approaching another simple limit of har- 
monic temporal shapes as illustrated by Fig. 9, where the amplitude is 9q = 0.5. 



This also follows directly from the linearization of equation ( 16 1 near zero 9 = 



4.2 Nonlinear case 

In this subsection, explicit analytical solution of the EP oscillator is obtained 
via transition to the action-angle variables, {9,9'} — > {/,</>}, of the generating 
model, k = 0. Re-scaling the time, p = eflGt, brings the effective Hamiltonian 



of EP oscillator (17) to the form 



H,, 



1 

2< 



?' 2 + tan 2 9 + fi cos 29) 



(23) 



where 9' = d9/dp is interpreted as a linear momentum, fj, = k 2 /(8G 2 ) is another 
parameter associated with the nonlinearity of Duffing oscillators k defined by 



(12), and the quantity G given by (15) must be treated as a fixed number 



Based on the exact solution (18), the action-angle variables are introduced 
as follows [221, 121 



arcsm 



V2I + 1 2 
1 + 



sin ( 



(1 + /) V2I + I 2 cos, 
v /l + (2/ + / 2 )cos 2 < 



(24) 
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Now canonical transformation ( 24 ) brings Hamiltonian ( 23 1 and the corre- 
sponding differential equations of motion to the form 



and 



dl 

dp 
d<p 
dp 



dHp 

' d(j> 



jit 1 + 1(2 + I) cos 20 

2 (TTTp 



1(2 + 1) . _ , 
MTi — ! — 7T9 sln 2</> 



(25) 



(1 + ^) 2 



0J 



2// 



(26) 



The essential advantage of new system ( 26 1 is that it becomes linear as \i = 



averaging can be also implemented as asymptotic integration of system ( 26 1 by 
means of the coordinate transformation {/, </>} 



while the original oscillator (231 still remains strongly nonlinear. The idea of 

ntegrat 
{J,V}: 

J(2 + J) 



I = J - fi 



ip - ft 



2(1 + J) 3 



cos 2tp + OO ) 



(J 2 + 2 J - 2) 
4(1 + J) 4 



sin 2^ + (a* 2 ) 



(27) 



Transformation ( 27 ) is obtained from the condition eliminating the fast phase 



tp from the terms of order fj, on the right-hand side in such a way that the new 
system takes the form 



dJ 
dp 

dp 



0(p 2 ) 
1 + J 



(1 + -0 3 



o( M 2 



System ( 28 1 is easily integrated as follows 



J = Jn 



(28) 



(29) 



1 + J 



(1 + J)3 



where Jq and ijj are constants of integration. 



Finally, substituting (29) in (27 1 and then (27) in (24), gives solution 6(p) 



whose effectiveness is illustrated by Figs. 10 and 11 at different nonlinearity 
and intensity of energy exchange levels, k and Pq, respectively. In particular, 
the transition from harmonic to sawtooth temporal shape of the phase 9 is due 
to the increase of the amplitude of energy partition index, Pq. The nonlinearity 
parameter n has certain effect on the curvature of lines as seen from Fig. 11. 
Note that, although the above procedure of asymptotic integration assumes the 
parameter h(k) to be small as compared to unity, the solution appears to be 
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very effective also under quite significant magnitudes of fi provided that the 
corresponding trajectory holds its symmetry on the phase plane 6-9'. However, 
the error of solution is growing as its trajectory becomes close to the separatrix 
of EP oscillator, which occurs under the condition n > 1; see Section [5] for 
physical interpretation. Such situation appears to be quite common due to 
highly sensitive dynamics near separatrix loop. 

Besides, expressions (27) and (29) point to the fact that, as the parameter 
fi increases, the monotonic growth of angle coordinate <j> can be broken. This is 
confirmed also by the direct numerical integration of equations ( 26 ) at different 
magnitudes of the parameter fi as illustrated by Fig. 12. The initial angle is 
<f> = tt/2 that associates with the extremum of the energy partition index P. 
The new action variable is fixed as J = 1.0, while the parameter fi is incremen- 
tally increased. As follows from the notations, this is eventually equivalent to 
incremental change of both the energy parameter k and the amplitude of en- 
ergy partition index, P . Detailed calculations show that above approximately 
fi = 3.559 (k = 1.4549 and Pq = —0.927062) the angle coordinate locks near its 
initial value <p = 7r/2, as a result trajectories on the planes represented by Fig. 
12 become closed. 



5 MODE LOCALIZATION 



The mode localization effect is shown schematically in Fig. 1. As mentioned in 
Section |3j the equilibrium points of oscillator ( 16 1 represent nonlinear normal 
modes of the original system. Although, in the case k^O, oscillator ( 16 ) is still 



solvable in quadratures, let us consider the following qubic approximation, by 
assuming that \6\ <C 7r/2, 



G - 



(8G 



2 ' k 2 ) 



= 



(30) 



As follows from ([7]), the equilibrium point 9 = of oscillator ( [30| c orresponds 
to the equal energy distribution, under which the original model |1| remains in 
one of its two symmetric nonlinear normal modes. So when the linear stiffness 



is positive, equation (30) has periodic solutions describing the energy exchange 
between oscillators ([5j near in-phase or out-of-phase mode. However, the equal 
energy distribution, associated with the equilibrium 9 = 0, becomes unstable 
if the linear stiffness is negative, G 2 — k 2 /4 < 0. In this case two new stable 
equilibria surrounded by separatrix loops occur near the unstable equilibrium. 
This indicates the onset of nonlinear local modes of the original system ([5| 
with a sustainable disbalance in the energy distribution despite of the perfect 
symmetry of system (JlJ. In terms of the present notations, the condition of 
negative linear stiffness can be represented in the form |21| 



nP 2 



2^/r 



p 2 

o 



2^1 



p2 
■M) 



= fl 



(31) 
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where Pq = - 
defined by ([£} 



sin 9{) and Qo — ~~ cos Aq the initial energy partition index as 



Condition (31) constitutes a necessary condition of localization because it 



does not guarantee that the dynamics is trapped inside one of the separatrix 
loops. The corresponding sufficient condition is obtained from the energy in- 



tegral of oscillator (30 1 in the form inequalities, first of which never holds for 
positive k, 



Qo > 



2 + 
2v/l 



^o 2 

-Pi 



9l 



(32) 



Qo < - 



2 - kP 2 



2^1 



pi 



.92 



Both estimates (311 and (32) are also valid locally, in the neighborhood of 
zero = 0, for strongly nonlinear oscillator (17). 



The above conditions (31) and (32) must be considered under the obvious 



constraint |Qq| = 1- Fig- 13 illustrates a relatively low nonlinearity case, when 
localization is impossible. The solid lines represent the boundary functions, 



introduced in (31) and (32), whereas the couple of dashed lines indicates the 



rectangular area within which solutions of inequalities (31) and (32) the above 
mentioned constraint. When the strength of nonlinearity n is increased, the 
line gi moves upward, whereas the line ji moves downward. When passing one 
through another at about k = 1, two small areas of localization occur as shown 
in Fig. 14. Note that, in both localization areas, the initial phase angle A lays 
in the neighborhood of zero. Therefore, according to Q, the localized modes 
branch out of the out-of-phase modes as the nonlinearity becomes sufficiently 
strong. 



6 THE BJJ FORMALISM 



6.1 Classical resonance of coupled oscillators 

The advantage of using the phase angle 9 develops through the form of oscilla- 



tor (16), which, in particular, admits clear qualitative analyses and allows for 
explicit analytical solution in terms of elementary functions when n = 0. This 



fact is essentially employed in Section 4.2 In this section, however, the energy 
partition index P will be used instead of the phase angle 9 in order to track 
the analogies with conventional approaches of macroscopic quantum dynamics 
as discussed in Section 6.2 below. Taking into account ^ and eliminating the 
ang le 9 from ^ and (jlijl, gives the coordinate transformation and the resultant 
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differential equations in the form, respectively, 



u 2 



V2 



1 

2 

-n 



= n 



K{l + P) cos 5 

-K{l + P) sin 5 
-K{l-P) cos(5 + A) 
^K(l-P) sin(£ + A) 



(33) 



and 



P 
A 



-efl^l-P 2 sinA= - 

w«p pcosA 



an, 



<9A 

_ dHeff 
dP 



(34) 



where 



H, 



eff 



-eCI 



en G 



kP 1 



sA\/l-P 2 



(35) 



and G is the result of substitution (|9| in (15). 

System ( 34 1 appears to have the form of B J J describing the interaction of two 



Bose-Einstein condensates in macroscopic quantum dynamics |27j |25] . Further 
details of this analogy are discussed below in the present section. Equations 



(34) are solved independently on the fast phase <5, which is determined then by 
integration from 



i 



k(1 + P) + cos A 



1 - P 
1 + P 



(36) 



Note that the effective Hamiltonian (35) can be obtained directly by sub- 
stituting (33) in ([3]) and averaging the result with respect to the phase S as 
follows 

4 _ 1 



Heff - 



2fl 



-eQn 



(37) 



Comparing the left-hand sides of equations ( 35 ) and ( 37 ) reveals the na- 



ture of integral (15), which is the constant mean value of the Hamiltonian, 
(H) s = const. It was already mentioned that, from the mathematical stand- 



point, the Hamiltonian system (34) and (35) with its simplifications is similar to 



that usually appears in macroscopic quantum dynamics as Josephson junction 
equations; see Subsection |6.2| for references and details. It known that such 
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systems are exactly solvable in terms of elliptic functions [27J , [55] . Eliminating 
the angle A in (35) by means of the first equation system (34), gives 



p2 



1 - P A 



2 



Ho 



(38) 



where H = -H e ff/(eCl)\ t=0 . 

Further manipulations leading to the elliptic functions are described in 



In particular, equation (38) gives the following quadrature 



P(0) 



dz 



(t) V( a2 + ^)( c2 -^) 



(39) 



where the parameters a and c are given by 

4 , 



2 2 

a c 



2 2 

a — c = 



12 \ l ~ H 0, 



(1 - kH ) 



Below, a series of diagrams (Figs. 15 and 16) provide qualitative descriptions 
of the dynamics in both A — P and Q — P planes. 



The evolution of topological structure of Hamiltonian ( 35 ) with gradually 



increasing nonlinearity levels is illustrated in Fig. 15. Similar diagrams were 
obtained in [25] for the case of two weakly coupled Bose-Einstein condensates. 
Classical nonlinear oscillatory chains were analyzed in somewhat different coor- 
dinates in |15j , and when investigating the intensity of energy exchange between 
parts of periodic nonlinear Frenkel-Kontorova and Klein-Gordon lattices based 
on the concept of limiting phase trajectories [2*5] . 

As mentioned at the end of Section [2} t he number of stationary points within 
the area, corresponding to rectangular ( 14 ), is increased by two when the nonlin- 



earity parameter goes above the bifurcation level k = 1. The new points, whose 
images are denoted by L in Fig. 15(c), correspond to the so-called local modes 
of the coupled Duffing oscillators branching out of the out-of-phase mode O; see 
also Fig. 1. As follows from the diagrams Fig. 15 (c-d), near the points L, the 
energy partition index P preserves its signature, in other words, the energy is 
accumulated in mostly one of the two oscillators despite of the perfect symmetry 
of the mechanical model. If the initial state (A(0),P(0)) lays outside but still 
close to the separatrix loop, which goes through the the image of out-of-phase 
mode 0, the dynamics is not localized any more but combines the effects of two 
stable local (L) and one unstable out-of-phase (O) modes. Such a combination 
becomes impossible however as further increase of the nonlinearity level k leads 
to the structural transition in the phase portrait such that the separatrix loop 
closes around the image of in-phase mode (I) rather than local modes (L). As 
a result, possible dynamics develop either near the in-phase mode or near one 
of the two local modes. Note that at high nonlinearity levels, n » 1, the 



contour lines of effective Hamiltonian ( 35 1 resemble the phase portrait of classic 
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pendulum as follows, for instance, from Fig. 15 (d). Based on such similarity, 
it is possible to simplify the effective Hamiltonian ( 35 ) as follows 



H efJ — > Hj = -eQ ( -kP 2 + cos A ) (40) 



The effect of such simplification on the contour lines, however, would lead to 
disappearance of stationary points L corresponding to the local modes. Hamil- 



tonian (401 gives the equation of pendulum whose angle is measured from the 



inverted (unstable) equilibrium 

A = (efl) 2 K sin A (41) 

Note that the conventional phase difference used in the macroscopic quan- 
tum dynamics is 4> = tt — A. The rigid pendulum analogy represents a typical 
reduction of SJJ equations describing the interaction of two weakly coupled 
superconductors, where $ is a relative phase, and P is a fractional popula- 
tion imbalance [25] . However, the BJJ in a double- well trap still requires the 
complete form of equations. 

Finally, using definitions ^ and (10), brings the diagrams of Fig. 15 to the 



form represented by Fig. 16. This gives a direct interpretation of the dynam- 
ics in terms of the energy partition and coherency indexes. In particular, the 
diagrams of Fig. 16 fully comply with the understanding of nonlinear normal 
modes as coherent motions of the system particles [30]. Namely, the normal 
modes correspond to the stationary points located on the vertical lines Q = ±1. 
In other words, when the energy partition index P is fixed, the oscillators co- 
herently vibrate either in-phase (Q = +1) or out-of-phase (Q = —1), according 
to the definition of normal modes. On the configuration plane u\U2, such type 
of motions is represented by pieces of lines passed by the system twice per one 
period. However, when the initial combination P(0) and Q(0) does not coincide 
with a stationary point then, during the vibrating process, both indexes slowly 
move along some trajectory on the plane PQ as shown in Fig. 16. 

6.2 Macroscopic quantum dynamics 

For comparison reason, let us reproduce few basic steps of the BJJ derivation in 
macroscopic quantum dynamics. Consider a two-level quantum system whose 
dynamics is described by the discrete version of nonlinear Schrodinger (Gross- 
Pitaevskii) equation [5] [23| (ft=l) 

dt d^* 

where tpj — ipjft) (j — 1, 2) are complex amplitudes of the wave function ip, 
and H is the system energy given by 

H = \{^ 2 + ^) + f (^i - Ml) - 1(^1 - Ml? (43) 
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Here V is the strength of coupling between two modes, while c is the strength 
of interaction between atoms, R characterizes the energy difference between two 
wells (bias). For the purpose of present work, it can be assumed that R = 



[3]; as a result, equation (42 1 gives 



dt 

.#2 

1 dt 



t/; 2 



1>l 



^2^*2)^1 
^2^*2)^2 



(44) 



In the literature, this model can be found in different physical contents and 
somewhat different notations [TS], [57], [25], [31], [10] • In general terms, it 
describes the tunneling effect in two coupled Bose-Einstcin condensates. 



The following complexification of transformation ( 33 ) is applied now 



^1 



ip 2 



-K(l 
2 V 



P) exp(Si) 



K(l-P) exp[(S + A)i] 



(45) 



\ip 2 \ 2 = 1, takes the 



where, the normalization condition for probability, IV'il 2 
form K = 1, and physical interpretation of the energy partition index P may 
depend upon the problem formulation. 



Substituting (45) in (44), gives equation for the fast phase 8, 



6 =2 P 



1 - P 
1 + P 



cos A 



(46) 



and the following classical Hamiltonian system for the canonically conjugate 
variables P-A, 



OH, 



P = -Vy/l - P 2 sinA = - 

oA 

A = 1 1 A -PcosA _ dH e ff 

VT^pt OP 

where H e ff is so-called effective Josephson Hamiltonian 

1 



eff 



H eff - - 



-cP 2 + V cos A Vl-P 2 



(47) 



(48) 



More often, however, the term Josephson Hamiltonian is used for a simplified 



version of (48), associated with the classical rigid pendulum as discussed at the 



end of Section 6.1 



Note that substitution of (|45|) directly into the original 

H,. 



Hamiltonian gives 

Ieff = 2H (49) 

Although the numerical factor in this relationship can be eliminated by re- 
scaling the time variable, nevertheless the presence of such factor moves trans- 
formation ( [45] ) out of the class of canonical transformations. This problem can 
be fixed [llj. however, at cost of losing the traditional form Josephson's system. 
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Despite the fact that the original system is quantum by formulation, the 



effective model ( 47 ) through ( 48 ) is usually qualified as a classical Hamiltonian 



system since the variables P-A can be simultaneously determined provided that 
it is true for their initial values. System (47) is equivalent to ( |34| under the 
conditions eSl = V and k = c/V. Substituting P(t) — — sm6(t) in (47) and 



(48), and eliminating the phase A, gives the effective conservative oscillator, 



similar to that was introduced for the case of coupled Duffing oscillators, 



H, 



c \ 2 tan ( 



eff 



cos 2 9 



1 9 

- -c 2 sin 29 = 



(50) 



where the number H e ff must be calculated by fixing a trajectory on the plane 
P — A or 9 — A before the oscillator is ( 50 ) is solved. 



Introducing the new time scale p = (H e ff + c/2) t and the parameter /i = 
(1/2) (1 + 2H ef f/c) "^enables one of using the solution in terms of elementary 
functions obtained in Section [472] for the case of coupled Duffing oscillators. 



6.3 Further analogies 

In terms of the macroscopic phase difference $ = <j) r — c/) l , corresponding to 
the right (r) and left (1) wells, and the atom number difference k — (Ni — 
N r )/2, Josephson type equations for the case symmetric double-well potential 
are obtained from Gross-Pitaevskii theory [57], 



$ = E c k 



Ej ^/l-Ak 2 /N 2 sin $ 
4k IN 2 



y/1 - 4k 2 /N< 



Ej cos $ 



(51) 



where N = iVj + N r is the total number of atoms; Ej and Eq and are called 
tunneling and charging energies, respectively, measured in units h. The charging 
energy arises from interatomic interactions, the tunneling energy determines the 
maximum current Ij = Ej. 
Introducing notations, 



2k = N t - N, 
N N t 

A = 7T - $ 



N r 



(52) 



brings system ( 51 ) exactly to the Hamiltonian form ( 34 )-( 35 ), after the following 



parameter substitutions 



ef2 



2Ej 
N 

N^Ec 

4 Ej 



(53) 
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Moreover, re-scaling the time variable as eClt = t in (34) and 2(Ej/N)t 
in (51 1, brings both systems to the same single parameter form [27J, 

dP 



dt 
d<f> 

~dl 



= -VI - P 2 sin$ 



= nP 



P 



V^p 2 



COS $ 



(54) 



where the strength of nonlinearity of Duffing system (12) k becomes equivalent 



to the parameter A [27] , the energy partition index between the Duffing oscilla- 
tors, P, is equivalent to the fractional relative population, and the phase shift 
between the Duffing oscillators, $, is the quantum phase difference between the 
right and left components. 

In order to use the results sections [3j |4.1[ and |4.2[ the fractional relative 
population P must be expressed through the phase angle 9. This leads to the 
integral 



He = -9 +-[ ff, jj 



NEc ^ t^ 2 9+( I ^\\os28 (55) 



1 



\ 8 



1 



where H e ff is the effective Hamiltonian of system (51) calculated on a fixed 
trajectory. 



Introducing the new time variable p and the action-angle variables (24) 



brings oscillator to its Hamiltonian form ( 26 ) , where 



(56) 



U 1 + ^n 

2 V NE C 



Therefore, analytical solution (24) through (29) obtained in Section 4.2 for 
the angle 9, gives 



P 



Ni~N r 
N t +N r 



sin 9 



V2I + 1 2 . 

sm < 

1 + / 



(57) 



where the action-angle variables, I and (f>, are given by (27) and (29) 



6.4 Second quantization 



As discussed at the end of Section 6.1 system (51) becomes equivalent to the 



classical pendulum under the conditions P < 1 and n >> 1 whose physical 



meaning is revealed by (53). As result, Hamiltonian of system (51) is reduced 



to that of the classical pendulum 



Hj = 



^k 2 - Ej cos$ 



k = - 



dHj 



$ = 



<9$ 
dHj 
~dk~ 



= -Ersm$ 



(58) 



= E c k 
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Except some technical advantages, such simplification does not seam to be 



crucial from the standpoint of classical dynamics since both systems (51) and 
(58) are integrable. However, in contrast to (51), the Hamiltonian of system 



( 58 ) has the quadratic form with respect to the 'linear momentum' k, which is 



important for the idea of second quantization 



Hj 



E c d 2 
' 2 <9$ 2 



E j cos $ 



(59) 



Further discussion and references on the transition from ( 58 ) to ( 59 ) can be 
found in [3] . Note that the effective Hamiltonian of energy partition oscillator 
( 55 ) is quadratic with respect to the 'linear momentum' 9 so that no further 



simplification would be required from the standpoint of second quantization. 



7 CONCLUDING REMARKS 

This work points to interdisciplinary links between the resonance phenomena in 
classical nonlinear dynamics and tunneling effects considered in the macroscopic 
quantum dynamics. It is shown that the classical dynamics of 1:1 resonance in- 
teraction between two identical linearly coupled Duffing oscillators is equivalent 
to the symmetric (non-biased) case of 'macroscopic' quantum dynamics of two 
weakly coupled Bose-Einstein condensates. The analogy develops through the 
BJJ equations, however, reduced to a single conservative EP oscillator impos- 
ing no additional assumptions. The derived oscillator is solvable in quadratures, 
furthermore it admits asymptotic solution in terms of elementary functions af- 
ter transition to the action-angle variables. The energy partition and coherency 
indexes are introduced to provide a complete characterization of the system 
dynamic states through the state variables of the EP oscillator. In particular, 
nonlinear normal and local mode dynamics of the original system associate with 
equilibrium points of the oscillator. Additional equilibrium points - the local 
modes - may occur on high energy level as a result of the symmetry breaking 
bifurcation. Further, it is shown that vibrations near nonlinear local modes 
of the coupled Duffing's oscillators are equivalent to the macroscopic quantum 
self-trapping (MQST) effect. As discussed earlier [35], there are two types of 
MQST, such as running-phase MQST and 7r-phase MQST. While BJJ captures 
both types of MQST, the SJJ rigid pendulum approximation ignores the 7T- 
phase MQST. Finally, it is noticed that, since the Hamiltonian of EP oscillator 
is always quadratic with respect its linear momentum, the idea of second quanti- 
zation can be explored in terms of BJJ without usual transition to the pendulum 



approximation. The quantum version of oscillator (55) allows for asymptotic 
limits with square well and sine-wave potentials, corresponding to the classic 
vibro-impact and harmonic oscillators, respectively. 
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APPENDIX 



As mentioned in the text, substituting ([6]) in ([5]) and then solving the set of 
equations with respect to the derivatives, gives 



K 



sKCL cos 9 sin(2(5 + A) + 



n 



fi cos 



sin 6-/2 cos 



sin(<5 + A) 



= e£l [sin A - sin 9 sin(26 + A)] - 
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Kfl 



/icos 



sin S + fi cos 



A = - 



2 s 



cos 8 cos(<5 + A) tan ( 



sec 9 



fi cos ( - — — I cos 5 + f 2 cos ( - + — ] cos(<5 + A) 



sin(S + A) 



(Al) 



O 



1 + e cos S cos(<5 + A) tan 



Kfl 



/1 cos 5 sec 



V 7T 

2 + 4 



where fi = U'{ui) and the coordinates itj (i — 1,2) are given by (JsJ) . 

System ( |60[ ) is still an exact equivalent to system ^ and represents a stan- 
dard dynamic system with a single fast phase, S. Applying the standard av- 
eraging procedure over the fast phase 6 on the right-hand side of system (60), 
gives equations (11). 
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" j£ Energy (Nonlineanty) 

FIG.l. Nonlinear normal and local modes of the coupled Duffing oscillators with 
symmetry breaking bifurcation generating two stable local modes from the out-of- 
phase mode. 
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FIG. 2. Geometrical interpretation of the coherency index on configuration plane 
in case of the energy equipartition, P = 0. 
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FIG. 3. Time histories of the energy partition index on the supercritical nonlinear- 
ity (energy) level K = 1.2 at different initial configurations of the model; see Fig. 4 
for the corresponding trajectories on configuration plane. 
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FIG. 4. Configuration plane trajectories on the supercritical energy level, K = 1.2, 
at different initial configurations of the model: (a) and (b) - the dynamics in small 
and large neighborhoods of the in-phase mode, respectively, (c) - the global dynamics 
around the unstable out-phase and both stable local modes, and (d) - the neighborhood 
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of one local mode; the corresponding energy partition index is shown in Fig. 3 (a) 
through (d), respectively, see also Fig. 5. 
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FIG. 5. Energy partition and coherency indexes on the supercritical energy level, 
K = 1.2, at different initial configurations of the model; see Fig. 4 for the correspond- 
ing trajectories on configuration plane. 
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FIG. 6. Configuration plane trajectories on the supercritical energy level, K = 2.0: 
(a) and (b) - the dynamics in small and large neighborhoods of the in-phase mode, 
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respectively, (c) - the neighborhood of local mode, and (d) - the drift over unstable 
out-of-phase mode and two stable local modes; see Fig. 7 for the energy partition and 
coherency indexes. 
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FIG. 7. Energy partition and coherency indexes on the supercritical energy level, 
k = 2.0, corresponding the dynamics illustrated by Fig. 6. 
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FIG. 8. Exact solutions for the beat dynamics of two identical linearly coupled har- 
monic oscillators and associated phase variables of the EP oscillator; highly intensive 
energy exchage. 
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respectively, for K = 0.591728 arid Pq = —0.430443; (a) phase trajectory, and (b) 
time history; p = eflGt. 
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FIG. 11. Analytical and numerical solutions represented by solid and dashed lines, 
respectively, for K = 1.41012 and Pq = —0.986701; (a) phase trajectory, and (b) 
time history; p = eflGt. 
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(a) 




FIG. 12. Angle trapping of the action-angle variables at high action levels (a), and 
its interpretation through the energy partition index (b); J = 1.0; [l = 0.5, 1.0,. ..,5.0. 
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FIG. 14. The plane of initial coherence versus energy partition showing two local- 
ization areas on higher nonlinearity level n. 
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FIG. 15. Contour lines of the effective Hamiltonian of the energy exchange equa- 
tions between Duffing oscillators at different nonlinearity levels: (a) quasi linear case, 
(b) post- bifurcation case k = 1.2 with center-saddle bifurcation near the out-of-phase 
mode, (c) strongly nonlinear multiple (>2) modes effect, and (d) strongly nonlinear 
effect of separatrix topological transition. 
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FIG. 16. Energy partition versus coherency index diagrams at different nonlinearity 
levels: (a) quasi linear case, (b) post-bifurcation case k = 1.2 with center-saddle 
bifurcation near the out-of-phase mode, (c) strongly nonlinear multiple (>2) modes 
effect, and (d) strongly nonlinear effect of separatrix topological transition. 
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